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Abstract 

The cross section data for ir° inclusive production in pp collisions is consid- 
ered in a rather broad kinematic region in energy y/s, Feynman variable xp 
and transverse momentum p?. The analysis of these data is done in the per- 
turbative QCD framework at the next-to-leading order. We find that they 
cannot be correctly described in the entire kinematic domain and this leads 
us to conclude that the single-spin asymmetry, A N for this process, observed 
several years ago at FNAL by the experiment E704 and the recent result 
obtained at BNL-RHIC by STAR, are two different phenomena. 
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High-energy single-spin asymmetries (SSA), usually denoted by A N , in 
pion inclusive production pp^ — > irX, have been observed for the first time 
more than ten years ago, by the Fermilab E704 experiment PjJ |2] at pi a b = 
200GeV/c (-^/s = 19.4GeV), in the beam fragmentation region. For clarity, 
we recall that A N is defined as 

_ do-(t) - dajj) 
AN ~ dtr(T) + da( I) ' (1) 

where (icr(T),(dcr(j)) is the n production cross section with one proton beam 
transversely polarized in the up (down) direction, with respect to the normal 
to the scattering plane. A striking dependence of An in Feynman variable xf 
was observed, as shown in Fig. ^for the neutral pion case, pp^ — > tt°X. Very 
recently, the STAR Collaboration |3| at BNL-RHIC has released new data for 
the same SSA at yfs = 200GeV, also depicted in Fig. ^ It is remarkable that 
these two sets of data cover approximately the same xp range and, assuming 
that they have the same dynamical origin, it might be tempting to conclude 
that this SSA is energy independent, to a first approximation. 

On the theoretical side, according to naive parton model arguments one 
expects An = 0, but several possible mechanisms have been proposed re- 
cently to generate a non-zero An- They are based on the introduction of a 
transverse momentum (kx) dependence of either the distribution functions, 
for the Sivers effect jlj or of the fragmentation function, for the Collins effect 
[Hj. These leading- twist QCD mechanisms have been used for a phenomeno- 
logical study of this SSA [6 j and higher-twist effects have been also considered 
[ZllHlin]- However a simultaneous description of A n and the unpolarized cross 
section has been ignored in all these works, a relevant point which has been 
already emphasized PHI and which will be crucial in the remaining of this 
paper. 

Our starting point is the analysis of the cross section data in a broad 
kinematic region in the framework of perturbative QCD (pQCD), at next- 
to- leading order (NLO). It is a common belief that it should be best working 
in the central region, namely, for a 7r° produced near 90°, in the center-of- 
mass system, although some authors have argued that the phenomenological 
inclusion of intrinsic kx effects leads to a better agreement between theory 
and data ^1] . The inclusive invariant cross section for the reaction pp — > tc°X 
reads 
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E T da/d 3 p n = J2J dx adx b f a/p (x a ,Q 2 ,fi F ) x 

abc 

fb/ P (%b, Q 2 , Vf) — /C ^ C ' - ' ^ da/ di(ab -> cX) , (2) 

where the sum is over all the contributing partonic channels ab — *> cX and 
da/dt is the associated partonic cross section. In our calculations the parton 
distribution functions (PDF) f a / p , fb/ p are the BBS set constructed in Ref. 
[T2"j and D n o/ c is the pion fragmentation function BKK of Ref. 2 ■ and 
fi' F denote two factorization scales and there is also a renormalization scale 
[Mr, associated to the running strong coupling constant a s . Our calculations 
are done up to the NLO corrections, using the numerical code INCNLL of 
Ref. [T5] . with fj,p = (jl' f = [ir = H- 

First we look at the cross section data at 90° for various energies, includ- 
ing data from ISR, Fermilab and BNL-RHIC, between y^i = 19.4GeV and 
200GeV, as a function of pt and the comparison with our calculations are 
displayed in Fig. El As it is well known, the agreement between theory and 
experiment depends on the choice of \i and this was extensively discussed, 
for example in Ref. [22]- For comparison we show the results for two scales 
li — pt and /i = pr/2 and we see that systematically the cross section for 
H = pr lies under that for \i = pt/2. In our case the best agreement occurs 
with n = pt for 200GeV and /i = Pt/2 for 52.8GeV and below. We get a 
beautiful description of the PHENIX data but we have some disagreement 
at lower energies; we underestimate the data by a factor two on average 
and more for E706, in particular, at low pt values. It is important to recall 
that the E704 Collaboration has also measured in this kinematic region 
6 = 90° and found = [IB], up to pt = 4.5GeV/c. The same trend has 
been observed recently by the PROZA-M experiment at pi a b = 70GeV/c [T7] 
and all these results are consistent with an earlier theoretical bound [THj . 

It is also crucial to note that in Fig. [TJ whereas the STAR data points 
are at a fixed angle 9 = 2.6°, this is not the case for the E704 data, for 
which 9° < 9 < 67°, due to the simple relation xf = 2pr/y / stan6 l and the 
fact that < pt > lies between 0.7 and 1 GeV/c. In this case, one can check 
that An = for the three lowest xp points which correspond to 9 > 15°, 

2 We have checked that the KKP fragmentation functions of Ref. ^1] give close numer- 
ical results 
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whereas is not zero only for the four highest Xp points which correspond 
to 9 < 15°. 

This observation led us to examine the unpolarized cross sections away 
from 9 = 90°, as a function of Xp, which combines the effects of pp and 9. We 
have considered the ISR data [20] which have the largest angle coverage, they 
are shown in Figs. El and 01 with our numerical calculations. It is striking 
to observe that the NLO pQCD results give a smaller and smaller fraction 
of the measured cross section, when going to smaller and smaller scattering 
angles. At yfs = 23.3GeV, which is very close to the E704 energy, for 9 = 15° 
the ratio Data/Theory can be one order of magnitude, even for // = pp/2, 
which means that another mechanism is at work and dominates, very likely 
related to soft processes of some specific nature. At -^i = 52.8GeV, we get a 
perfect description of the 9 = 53° data, for fi = p T /2, but again the relative 
size of the soft contributions increases in the forward direction, where the 
disagreement between data and pQCD gets larger. The failure of the NLO 
pQCD to reproduce the ISR small angle data is certainly not due to the 
breakdown of the convergence of the pQCD series in this region. Indeed 
the " K factor ", defined as usual as daNLo/d&LO, where LO stands for 
leading-order, is larger for small angles. We have also observed that there 
is a considerable reduction of the scale dependence, going from LO to NLO 
which shows an improvement in the perturbative stability. Therefore next- 
to-next-to-leading or higher pQCD corrections are not expected to account 
for a factor ten or so in the ratio Data/Theory. 

Finally, we show in Fig. the NLO pQCD calculations near the forward 
direction at ^ = 200GeV and one notices that at this energy, the 9 = 2.6° 
STAR data is very well reproduced by the solid curve (fi = pp). For the sake 
of completeness we also give our prediction (dashed curve) at 9 = 4.2°, com- 
pared to the preliminary data released recently by the STAR Collaboration 
jHj- Now let us go back to the SSA measured by E704 and STAR shown in 
Fig. CI In both cases, the data lie near the forward direction, but in the 
case of E704 it cannot be attributed to a pure pQCD mechanism, because 
the measured unpolarized cross section can be one order of magnitude higher 
than the pQCD result (see Fig. EJ). On the other hand, it is very likely to be 
so in the case of STAR which has measured the cross section in full agree- 
ment with pQCD. Therefore this analysis strongly suggests that the two sets 
of SSA data, we have considered here, are the manifestation of two different 
phenomena, in contrast with the implications from Refs. |d1 I7| I23j. 

To conclude, if pQCD does not generate a SSA in the 9 = 90° region, as 
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shown by E704, we dare to predict that PHENIX should also find = 0, a 
result, indeed perfectly consistent with the recent rigorous positivity bound 
[21], which implies \A^\ < 1/2. It will be very interesting to know what the 
data analysis will give, when completed in the near future 3 . 

Note added : After this paper was completed, L. Bland and A. Ogawa have 
compared PYTHIA calculations of the invariant cross section with STAR and 
ISR data and they confirm our conclusions. 

Acknowledgments: We thank W. Vogelsang for checking our numerical 
results and a very useful discussion. We are also grateful to M. Begel, G. 
Bunce, L. Bland and N. Saito for interesting comments and suggestions. 



3 This has been confirmed by very recent preliminary data |25| . 
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Figure 1: The single-spin asymmetry An as a function of Xp, at two different 
energies. The data are from Refs. [21 E] . 
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Figure 2: Ed 3 a/d 3 p at 90° and various energies, as a function of px- Data 
are from Refs. ^TJ QHJ EH1 Elj and the curves are the corresponding NLO 
pQCD calculations with // = pt (solid lines) and \i = pr/2 (dotted- dashed 
lines). 
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Figure 3: Ed 3 a/d 3 p at = 23.3GeV, as a function of xp for two different 
scattering angles. The data are from Ref. [20j and the curves, solid 6= 
22° and dashed 9= 15°, are the corresponding NLO pQCD calculations with 
(jl = pt- The dotted-dashed curves are for n = p T /2. 
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Figure 4: Ed 3 a/d 3 p at y/s = 52.8GeV, as a function of xp for three different 
scattering angles. The data are from Ref. [201 and the curves are the corre- 
sponding NLO pQCD calculations with fi = p T . The dotted-dashed curves 
are for \i = px/2. 
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Figure 5: Ed 3 a/d 3 p at y/s = 200GeV, as a function of xp. The solid and 
dashed curves are the NLO pQCD calculations, with /i = p?, at two different 
angles and the data points are from Ref. The dotted-dashed curves are 
for \i = Pt/2- 
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